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1 Introduction 


The chiral constituent quark model is a mixture of the potential model and the 
cr-model which contains spontaneous chiral symmetry breaking in a dynamical way. 
The quarks get their masses via the interaction with the scalar field, which develops 
a vacuum expectation value. In a previous paper the thermodynamics of this 
model has been investigated, in particular the restoration of chiral symmetry with 
temperature. It is important to understand the possible evolution of this model in 
QCD in order to handle it correctly at finite temperatures. For high resolution, 
i.e. short distance processes fundamental QCD with quarks and gluons is the most 
efficient theory. A typical scale associated with perturbative QCD is A > 1 — 1.5 
GeV. RHIC and LHC physics for particles with p_L > A will be dominated by such 
processes. In nucleus nucleus collisions, however, scattering at smaller momentum 
scale will be non negligible. Since the strength of the QCD interactions increases 
with decreasing momentum transfer, characteristic qq bound states will form and in¬ 
fluence the dynamics at larger distances. A typical resolution where these processes 
start to become important is A^ 5 b = 0.5 GeV [|]. Below this scale the vacuum 
changes and acquires a quark condensate and/or meson condensate. The sponta¬ 
neous breaking of chiral symmetry will persist down to A ^qcd- Around this scale 
the confining gluon configurations make themselves felt via confining forces. In the 
interval [Aqcd, -^xSb] dynamics is governed by constituent quarks interacting 
via pions and sigma mesons. One should keep in mind that the effective couplings 
of these mesonic degrees of freedom to themselves and to the quarks are induced by 
the underlying physics of the quarks and gluons. When one immerses the quarks 
and gluons into a heat bath, it is to be expected that the induced couplings will 
be temperature dependent. In the following paper we calculate the temperature 
dependence of the couplings in the linear sigma model. We apply the Heat Kernel 
technique to integrate out the fluctuations of the fermion variables. 

We use the constituent quark model as an effective theory in the region between 
^QCD ~ 200 MeV and A^sb ~ 0.5 GeV with constituent quarks and confining glu¬ 
ons, which are neglected in this paper If the integration over quark and meson 
fields is done in one step, then one obtains a first order phase transition contrary 
to our expectations for the sigma model and QGD for zero quark masses. This 
calculation therefore presents an intermediate step to illustrate the technique of the 
heat kernel which will be finally applied in the context of the renormalization group 
equations to get a more accurate description of the second order phase transition. 
In the course of developing this technique we solve the problem of quark condensate 
in the constituent quark model. The natural order parameter in the linear sigma 
model is the vacuum expectation value of the scalar field. We show how we can 
relate the vacuum expectation value of the scalar field to the quark antiquark con¬ 
densate. In addition we calculate TTTr-scattering as a function of temperature. This 
calculation shows the slightly different physics concept underlying the constituent 
quark model and the Nambu-Jona-Lasinio (NJL) model in a physical observable. 
The TTTT-scattering amplitude can be tested via transport properties. We also show 
that no new renormalization constants or cut offs are necessary in the constituent 
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model contrary to the NJL model, where the meson loops generate additional sin¬ 
gularities which have to be regulated by new phenomenological cut off parameters. 
The heat kernel technique allows to restrict the effect of quarks of a hxed mass to 
the bosonic vertices. That means we propose to ignore the free fermion loop contri¬ 
bution of unconhned quarks in the free energy and also to freeze the quark mass in 
order not to allow artihcially long distance meson meson interactions, which would 
be cut off by the conhning forces between quarks and antiquarks below the critical 
temperature. In the constituent quark model it is assumed that the meson degrees 
of freedom persist beyond the chiral phase transition. Since the gluonic forces are 
especially strong in the tt and a channels these mesons may not dissolve but remain 
active degrees of freedom beyond the chiral transition. 

The outline of the paper is as follows: In the next section 2 we derive the heat 
kernel method at zero temperature and apply it to the constituent quark model at 
hnite temperature in section 3. We calculate the free energy, which is the necessary 
thermodynamical tool for the analysis of a phase transition and study the chiral 
transition. In section 4 we evaluate the quark condensate and the Trvr-scattering 
length as functions of temperature explicitly. At the end we conclude and compare 
our results to other calculations. 


2 The Heat Kernel Method for T = 0 


We consider the chiral constituent model with quarks, a and tt mesons. At zero 
temperature in Euclidean space the partition function or the generating functional 
without external sources A is given by 

Z[A = 0] = j VqVqVaVn exp{— J {Cp + Cb)} (1) 

with a fermionic Cp and a bosonic Cb part. 

Cp = q{x) {jpdp + g {a + iTnj 5 )) q{x) (2) 

+ -^(cT2 + ff2) + ^(a2 + 7f2)2 . (3) 

The Yukawa coupling of the constituent quarks to the mesons is denoted by g. 
The parameters of the linear a-model at T = 0 are chosen in the same way as 
in ref. |]I|. The mass squared /Xq is positive reflecting the spontaneously broken 
ground state. The minimum of the potential U{a) lies at (ao) = 0.093 GeV. The 
light constituent quarks have a mass of m = 300 MeV, the cr-mass mg- = 700MeV, 
Po = (0.495 GeV)^, A = 28.33 and g = 3.23. These couplings are effective couplings 
obtained after the integration over high momentum quark and gluon degrees of 
freedom. Therefore we must not reintegrate these high momentum modes at zero 
temperature, since their effect is already included in the couplings. In the following 
we demonstrate the heat kernel technique at zero temperature to set up the notation 
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we are going to use at finite temperature. Also in a later paper we are investigating 
the renormalization group flow with the same method [O . 


By means of a formal integration over the quarks one obtains a non-local deter¬ 
minant, which is in general complex 

Z[A = 0 ] = J VaVif det {'yEdE + g {c + ifTT'y^)) exp{— J cI'^xCb} 
det {'JeOe + g {o'+ iff'y^)) = J 'Dq'Dqexp{- J (I'^xCf}- 

This determinant can be expressed as a logarithm and in the following we consider 
the real part of the effective action (the modulus of the determinant) 


= —Tr In T)T)+ 


(4) 


with the elliptic operator 

DD+ = -dl + g'^MM^ +g^^{d^M^) 

and M{x) = a(x) +iff(x)'y 5 . The symbol ’Tr’ is dehned by a functional trace over 
an operator O 


TrO = j d!^xtr{x\0\x) 


and tr is reserved for the remaining trace over all inner spaces i.e. spin, positive and 
negative energy solutions associated with the gamma matrices, isospin and color [ffl . 


Since DD^ is a positive dehnite operator it is possible to rewrite the logarithm 
as an Schwinger proper time integral. Then we hnd for the real part of the effective 
action in d = 4 dimensions 




lim - 

A^oo 2 


dr _ 

/ — Tre 

Jl/bP- T 


— T 


DD+ 


1 dr 
hm - — 

A^oo 2 Jl/A2 r 


d'^xHo{x, x; t) 


where the second line dehnes the diagonal part of the zero temperature heat kernel 
Ho {x,x-,t). In general the proper time integration diverges for vanishing tempera¬ 
ture therefore we introduce a proper time regulator A. The order of magnitude for 
A is estimated in ref. @ from the meson spectrum and lies in the range of about 1 
GeV. Physically one can interpret A as a scale at which chiral symmetry breaking 
starts. 


The non-diagonal part of the heat kernel Hq{x, y; r) can further be evaluated by 
including a plane wave basis p. We get 


Ho{x,y,T) 


tre {x\y) 


tr 


tr 


^ ^ -iky -rDDt -ikx 

{2nY^ 


(5) 
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which yields the diagonal part {y x) of the heat kernel 


Ho{x,x]t) 


tr 


(27r)^ 


^-T{k^-2ikfj,d^+DDt^ 


1 . 


Note, the unit operator on the right is to be associated with other possible functions 
on which the derivatives with respect to x can act. The product rule of differentiation 
can be implemented by the substitution (9^ —+ in the operator DD^. After a 
subsequent rescaling of the momenta k^ kfJ,/^/T the following expansion in powers 
of the proper time r is found 


Ho{x,x;t) 


tr 

(47rr)‘^A 



( 6 ) 


The momentum integration can be done using the relation 


dk^k 

(i/2* 


TT' 


kfMl kf^2 ■ ■ ■ kf^2m 




Om E *^A‘<t( 1 )Mct( 2 ) *^Mct( 3 )M<t( 4) 

(7G5(2m) 


•■5 


M(T(2m—l)Mcr(27n) 


with (t(1) < (t(2) , (t( 3) < (j(4) , ■ ■ ■ and (j(l) < (j(3) < (t( 5) ■ ■ ■ . All in all there 
are (2m — 1)!! permutations. For the constituent quark model the corresponding 
expression for the heat kernel in d = 4 dimensions looks like 


Ho{x,x,t) 


tr 

(Anry 


^2 


with n := g'^MM^ and V := g'y ■ [dM^). A mass term md in the exponent has 
been introduced to guarantee the convergence of the proper time integration. It is 
chosen equal to the quark constituent mass in the vacuum. Quantum fluctuations 
of the theory for short distances enter in the r —0 region of the r integration while 
the behaviour of quantum fluctuations at long distances contribute in the limit r —> 
cx). Therefore through the introduction of a mass term m^ no infrared divergences 
for r —> CX) emerge. The mass term m^ is an IR regulator and is not uniquely 
determined. We checked the sensitivity of our later results with respect to variations 
in m^ and found that in a reasonable interval there is no strong dependence. 


The heat kernel method makes the fermion functional integration doable by 
expanding the exponential in the proper time variable up to the second order. Per¬ 
formance of the inner traces yields the following result for the effective action 

(•oo ftr 1 9 

= AN, lim / — -— 

a^ooA/a 2 r (47rr)2 

1 - 'r + y 

+g'^{{d^ay + {dfj,ny) + m'^ - 2m'^g‘^{a‘^ + + g^d‘^{a‘^ + ^^)} • 
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The proper time integration leads to an incomplete gamma functions in the fermionic 
action 


= 






E 

fc =0 


m 


2k 


Tr hi 


where Trhk symbolizes a short hand notation for the integrals resulting from dif¬ 
ferent powers in r containing the helds and their derivatives |Q. One observes that 
the proper time expansion leads to a gradient expansion, since the proper time has 
dimension mass~'^. The first three coefficients of this expansion are always divergent 
for A —cxo, therefore a regularization scheme should be applied at this stage. 


3 The Heat Kernel Method at finite Temperature 


Before we apply the heat kernel method to the constituent quark model we want to 
show how to extend it to finite temperature. 


The main difference in finite temperature field theory compared to zero temper¬ 
ature field theory is the change of the continuous time integration of the effective 
action to a hnite integration in Euclidean space from zero to /3 = 1/T. We can 
generalize the definition of the heat kernel at zero temperature cf. Eq. (^) to finite 
temperature 

{x\x') p 

= - x) (7) 

where '[3' stands for the effect of finite temperature in the Fourier representation 
of 5^{x — x'). This means one has to insert the correct boundary condition of the 
fermions (or bosons) into the Dirac delta function. For the quarks we have to take 
anti-periodic boundary conditions in Euclidean time, i.e. the Fourier representation 
has only half integer Matsubara frequencies. With the help of a generalized theta- 
transformation one can End a connection between the zero temperature and the 
hnite temperature heat kernel for fermions |P 


Ht{x,x,t) 


Ho{x,x,t) 


1 


OO 

n=l 


( 8 ) 


This result is remarkable, because the hnite temperature heat kernel separates into 
a sum of the zero temperature heat kernel and a temperature dependent piece, 
which again is a product of the zero temperature heat kernel times a temperature 
dependent function. 


Since we are interested in the pure temperature contribution of the ehective 
action we consider only the diherence between the hnite temperature and zero- 
temperature heat kernel. The temperature dependent ehective action is then given 
by 





Ho{x,x;t)} . 
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The vacuum divergences in both kernels cancel each other and no new divergences 
are generated at hnite temperatures. Therefore we do not need the regulator A 
anymore and take the proper time integration over all positive real r. With the 
result of Eq. (Q) we find 

1 r /*oo /~J 

= - d^x 

2J Jo T 


In the previous section the expansion of Ho{x,x;t) was shown and by using the 
formula 

u/2 


= 2 ( - | K,{2^) 


7. 


where K^{x) denotes modihed Bessel functions, we hnally obtain 
“ d^x 


^eSl = 


Nr. 8m^ „ / ^ X Nr 4m „ , ^ , r 9/ 9 ^9 n 9-, 

F2(^m) - ——Fi(^m){/(cT^ + tt^) - m^} 


Nr 


27r^ 


2ti‘^ y 2 tt'^ (3 

Fo{(3m) {m^ + + {d^f) + /(cx^ + Tf^)^ 

1 


-2m^g\ay + -g‘^d\ay 

o 


Here we define 


F.(M := E(-l)' 


, Ki{nf3ni) 


n=l 




for i = 0,1, 2 . 


(9) 


This result is substituted into the partition function Eq. (|l|), which then becomes 
only a functional of the mesonic degrees of freedom alone. 


Z= DaDn exp 


+ X (10) 


exp — / 


y{T) 


--AT){{d,ar + {dyf)}- 


A + y + 


A'(T) 


2^2 


+ A 


The fermionic integration has given a prefactor which comes from the quark loop at 
hnite temperature. This prefactor is due to the free motion of quarks, which we think 
is inadequately described by the linear sigma model without conhnement. Therefore 
we will not include it in the equation of state. We take, however, the fermionic effects 
on the mesonic couplings seriously. Due to the freezing of the infrared parameter m, 
the fermion loop contributions arise from short time huctuations in the thermal 
system, where quark conhnement is not relevant. After sorting one encounters 
the following temperature dependent wave function renormalization constant and 
couplings up to second order in the proper time expansion: 


z{T) 

AiT) 

A'(T) 


1 

Ao + + 

O AT 


8 Nr 8m 
(47r)2 13 


g^Fi (Pm) 


( 11 ) 

( 12 ) 

(13) 
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The first terms on the right hand sides in Eqns. ( pT| - |T^) are the T = 0 terms coming 
from the bosonic action Eq. (^. The finite temperatnre corrections are generated 
from the heat kernel. Since the heat kernel expansion is an one loop integration, 
its result corresponds to the results in perturbation theory which gives temperature 
dependent modifications of the couplings. The wave function renormalization comes 
from the graph where a quark radiates a meson. The mass of the mesons is changed 
by the meson polarization operator in quark antiquark pairs to order g'^. The meson 
meson scattering amplitude gets modihcations from the box graph in fourth order 
of the coupling. The numerical values of the couplings are shown in Fig. as a 
function of temperature. One sees that p' decreases and becomes zero at T 185 
MeV. The coupling A' and the wave function renormalization z also become smaller. 
Note the wave function renormalization z does not vanish at the chiral transition, 
which means the pion has not dissolved at this temperature. The probability to find 
a free pion at T 185 MeV, however, has decreased to 50%. 



Figure 1: couplings A', g.' and 2 ; as function of the temperature. 

We also apply the presented heat kernel method to evaluate the quark condensate 
{qq) at hnite temperature in the chiral limit and to calculate TTTr-scattering for finite 
pion masses. These calculations are possible by adding an additional source term 

2 

2Agg = A ^ 

i=l 

to the fermionic Lagrangian Eq. and then differentiating with respect to A. Thus 
the parameter A plays the role of a current quark mass and breaks chiral symmetry 
explicitly. We again calculate only the difference between the finite temperature 







heat kernel Ht and the zero temperature Hq. Therefore we have to include in our 
bosonic action at zero temperature an explicit symmetry breaking term, arising from 
the integration over the quarks with finite masses at zero temperature which we do 
not execute explicitly. This term will receive a temperature dependent modification 
in our calculation. In addition another term cubic in meson fields will arise at 
hnite temperature through the fermion integral. The symmetry breaking term in 
the linear sigma model at T = 0 is given by Cc = —ca with a value of c which 
reproduces the physical pion mass. 

ca = -mlf^a. 


Using the relation rri^f^ = —2A(gg)o we get 


c = A 



=: Ac 


Employing the heat kernel technique we obtain the final effective bosonic action at 
T 7 ^ 0 including explicit symmetry breaking due to finite quark masses A. The 
meson wave function renormalization constant has been scaled out. 


S'[cr, ^] = J I const + ^ + 

+ts{T)Aa ^ + ti(T)Acr| 

with additional temperature dependent functions entering the symmetry breaking 
terms. 


ti{T) 

tsiT) 

A"(T) 


-2(/r-(T)^-/ig) + 5 

z(T)V2 

4(A'(T) - Ao) 

z(T)3/2 

z{T) 

A^(T) 

z{Ty 


(14) 


There are two new temperature dependent loop corrections in this formula propor¬ 
tional to the external source A. The sigma tadpole term tiAa gets a temperature 
dependent fermion loop correction. A temperature dependent Yukawa vertex be¬ 
tween two mesons and the sigma meson is generated due to the quark loop. The 
term t^A adds to the Yukawa vertex produced after spontaneous symmetry break¬ 
ing, when a develops a vacuum expectation value. In Table 1 , we give the numerical 
values of /x". A", fiA/(/^m^) and t^A. Both terms depending on the explicit sym¬ 
metry breaking A influence the order of the phase transition. Previously 


we 


did not allow such temperature dependent terms in calculations of the order of the 
phase transition. Especially a third order term would add to the determinant term 
in Flavour-S'f/(3) from instantons which is nine times larger at Tc 
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T [MeV] 

0 

40. 

80. 

120. 

160. 

Tc 184.8 

/i" [MeV] 

495. 

495. 

493. 

469. 

362. 

0. 

A" [1] 

28.3 

28.4 

30.1 

36.9 

53.0 

72.9 

Ati/Am^ [1] 

-1. 

-1. 

-0.97 

-0.81 

-0.42 

0. 

Afs [MeV] 

0. 

-0.1 

-7.4 

-38.2 

-116.4 

-219. 


Table 1: Couplings in the effective Boson action S'[ct, vf] for different temperatures. 


4 Quark Condensate and TTTT-Scattering 


The next task is to track the behaviour of the vacuum expectation value of the a 
field. When we expand the effective action around the spontaneously broken ground 
state (To = const, we get up to a normalization factor J\f the partition function with 
the tree contribution to the effective action. This action contains all the effects of 
the temperature fluctuations of the fermions besides the free energy of free quarks 
and antiquarks. These effects suffice to produce a chiral phase transition. 

,/'2 \// 

S[ao\ = {(3V) (tiA cto - —al + + —ao) 

In mean field approximation we determine from this tree action the quark condensate 
in the chiral limit as a function of temperature. 


2{qq) 


MF — 


1 d\nZMF{A) 


(3V 


dA 


— + tsCTg. 


(15) 


A=0 


Minimizing the action with respect to a one finds that (Jq and the resulting quark 
condensate are smooth functions of temperature (cf. Fig. which vanish at the 
chiral transition temperature where vanishes. The phase transition is at a higher 
temperature compared to the calculation in ref. , because the effects of the quark 
fields enter only via the intermediate coupling of the meson fields. If one would set 
m = gao in the prefactor in Eq. (pTf) , a variation of this part would give a faster 
decrease of the vacuum expectation value of (Jq. But this procedure is questionable 
since the quarks are not free below the transition temperature. 


In order to estimate the mesonic effects on the quark condensate we include 
their thermal fluctuations in the partition functions. In the linear sigma model 
the integration over the bosonic degrees of freedom is much simpler than in the 
NJL model p. The zero temperature bose fluctuations are already included in 
the meson coupling constants at zero temperature. Therefore we use a simple one 
loop approximation for the boson integration. Including mesonic terms we obtain 
the following partition function, where the integrals are to be understood at finite 
temperature with the mean field masses. 

2 = _^^-SN-irrln(-a2 + Af|)-|rrln(-a2 + Af|) 
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Temperature [MeV] 


Figure 2: Normalized condensate {qq)MF and {qq) with bosonic fluctuations as 
function of the temperature. 


= -/i"^ + 3 A'Vq + htsAcTo 

+ A Vq + 2t3A(To (16) 

The value ag entering these equations is determined in mean held from the minimum 
of the bosonic tree action. Consequently, ag depends on the explicit symmetry break¬ 
ing parameter A. The quark condensate in this approximation contains fermionic 
and mesonic effects. If the conplings ti and tg wonld not depend on temperatnre, 
the qnark dynamics wonld be invisible in the variation of the qnark condensate with 
temperatnre. 

‘^{qq)T = ti(T)ag + t3(T)ag -|- /cr(o'o, T) + I-^i^ag, T) (17) 


where the bosonic integrals are defined as follows: 


(27r)3 


M,T dMl 

(9A 


E 

A=0 n=l 


K^inMjT) 


and 

h{(yo,T) 


3 I" (Pp 1 1 dMl 

2 7 (27r)3 + Ml - 1 '9A 

T2 dMl 

^ ^ A=o ’ 
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In the last equation we have to take into account the pion with zero mass. To specify 
the result in the chiral limit we list the masses and cxo held derivatives 


dM^ 


dA 

dM^ 


A=0 

A=0 


dao 


dA 


A=0 


+ {2X"ao) 


dMl 


dan 


dA 


A=0 


dA 


A=0 


tl + StsCTg 


Ml 


A=0 


(18) 


In Fig. 2 we show the resulting quark condensate including the bosonic thermal 
huctuations. It decreases faster and leads to a hrst order transition. We are aware 
that the meson one loop approximation is not sufficient to generate a second order 
phase transition for a pure mesonic (9(4)-theory. In the present paper we are mainly 
interested in the effect of quarks on the mesonic action, in a further paper we are 
going to improve our calculation by combining renormalization group methods with 
the heat kernel technique [|^. 


The TTTr s-wave scattering amplitude vanishes in the zero momentum limit for 
massless pions. Therefore we have to use the masses of a- and vr-mesons in the form 
with hnite symmetry breaking A ( cf. Eq. (pffiD). In addition we need the Yukawa 
coupling ^fCTTTTT from the bosonic action after symmetry breaking: 


9o 


— A'Vo + t^A 


(19) 


The s-wave projection of the Trvr-scattering amplitude in the isospin zero channel 
then has the following form 


roo(s = ml) 



Ml-AMI Ml\ 


( 20 ) 


Inserting the temperature dependent masses Eq.(|T^ and coupling constants Eq.(p)4|) 
we hnd a divergence of the s-wave isospin zero scattering length (cf. ref. [^]) shortly 
before the chiral phase transition. As long as one treats the bosons very naively in 
mean held or one loop approximation there comes a point where the a mass becomes 
degenerate with twice the pion mass, which leads to an exploding scattering length 


floo 


Ooo 


^00 

32ttM.^ 


( 21 ) 


In Fig. y we plot the result of our calculation near Tc and compare the result with 
a simple extrapolation of the Weinberg scattering length 


ooo 


IrriT, 

3^ 


to hnite temperature. 


A more careful investigation of pion scattering at hnite energies shows that the 
ehect of the singularity is diminished when one averages over the cm-momentum of 
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Figure 3: Scattering lengths ago as functions of the temperature. 


the pions, as the heat bath only contains hnite momentum pions. The total effect is 
almost completely washed out, if one includes an average over the unitarized s-wave 
amplitude, which contains TTTr-rescattering effects. It is therefore very important 
to include the ring sum or an equivalent damping mechanism to control the pion 
interactions near Tc. The scattering length is a good illustration for the inadequacy 
of the mean held or one loop approximation in the presence of Goldstone Bosons. 
It remains to be seen how far away from these effects are important in the or¬ 
der parameter, since in the physical mass case the second order phase transition is 
smeared out. Also a quantitative comparison of the ring or large Nf approximation 
with the renormalization group equations is interesting. The heat kernel technique 
allows without any additional assumption about the form of the infrared regulator 
to calculate the evolution equations. At lower temperatures the presented coupling 
constants can be included in the constituent quark model calculation with conhne- 
ment in ref. to take into account the effects of thermal quark huctuations on 
the couplings of the sigma model. The simplicity of NJL model estimates remains 
unchallenged by the linear sigma model. If one wants to tackle higher order effects, 
however, our calculation shows that the constituent quark model is by far superior 
in its ability to include meson dynamics. The physical question whether mesonic 
degrees of freedom are modeled adequately by this model near has to be answered 
by calculating susceptibilities not only in the scalar and pseudoscalar, but also in the 
vector channels. These results can be compared with numerical lattice calculations 
but also with stunning experimental data from dilepton spectra [O . 
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